Nonequilibrium Campbell length: probing the vortex pinning potential 
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The AC magnetic penetration depth X{T,H,j) was measured in presence of a macroscopic 
DC (Bean) supercurrent, j. In single crystal BSCCO below approximately 28 K, A(r, II,j) exhibits 
thermal hysteresis. The irreversibility arises from a shift of the vortex position within its pinning well 
as j changes. It is demonstrated that below a new irreversibility temperature, the nonequilibrium 
Campbell length depends upon the ratio j/jc- X{T,II,j) increases with j/jc as expected for a 
non-parabolic potential well whose curvature decreases with the displacement. Qualitatively similar 
results are observed in other high- Tc and conventional superconductors. 

PACS numbers: 74.25.Nf, 74.60.Ec, 74.60. Ge 



The AC penetration depth A is an extremely sensi- 
tive probe of the vortex state. Models for the small- 
amplitude AC response usually assume a uniform distri- 
bution of vortices and a parabolic effective potential well 
iQ-^. However, experiments often probe a nonequilib- 
rium flux profile. In this case, A may depend not only 
upon H and T but also upon the Bean supercurrent j 
arising from a gradient of vortex density. In this paper 
we show that a finite j produces a hysteretic AC re- 
sponse whenever the pinning potential is non-parabolic. 
In contrast to most experiments, the irreversibility dis- 
cussed here occurs with temperature while the applied 
field is held constant. Thermal hysteresis is observed in 
several materials, but we focus on BSCCO. In all our 
experiments Hac 5 mOe) <C Hdc ^ Hc2 where the 
excitation field does not modify the average vortex dis- 
tribution and we can ignore vortex core contribution to 
A. 

Vortices transmit the perturbation caused by a small 
AC field as either compressional or tilt waves, depending 
upon the geometry of the experiment |^ . If the ampli- 
tude of the AC field is small, the elastic restoring force is: 
F = —au, where u is the vortex displacement caused by 
a small AC current and a is the Labusch parameter. The 
latter incorporates, self - consistently, both the pinning 
and the elastic forces. The AC penetration depth is then 
given by A^ = Xj^+X^ where Ac is the Campbell penetra- 
tion depth . The Campbell length is Ac = \fCcxfoi, 
where Cxx is the appropriate elastic modulus [C\\ for 
compression or C44 for tilt, C44 « Cn w _B^/47r) 
Unlike earlier versions based on the local pinning force 
these formulae are also good approximations for a 
non-local elastic response with dispersive a(fc) It is 
also possible to incorporate fiux creep and flux flow effects 
. These effects, however, do not result in thermal ir- 



reversibility of X{T, H, j). In this paper, we explore the 
simplest generalization of the elastic model. In the pres- 
ence of a macroscopic DC Bean current, j, the Lorentz 
force biases the mean vortex position away from the pin- 
ning potential minimum. The Campbell depth is then 
determined by the small - amplitude response near this 
new bias point. If the dependence of the pinning poten- 
tial upon vortex displacement u is non-parabolic, then 
changes in either j or the critical current jc will result in 
changes of the Campbell length. If j relaxes via flux creep 
or jc changes with temperature, the Campbell length will 
then show a small, but observable thermal hysteresis. 

Let V{x) describe the shape of the pinning poten- 
tial with X = r/vp, where r is the vortex displacement 
and Tp is the effective pinning radius, rp ^ ^ for a 
quenched disorder pinning at low temperatures ( see Fig. 
|l|). The Lorentz force due to j causes a dimensionless 
vortex displacement xq determined from the condition: 
dV/dx — Bj/rpC — ao(j/ic), where ao = Bjc/cvp. Upon 
application of a small AC field, the restoring force acting 
on a vortex in the vicinity of xq is F{u) ~ ~au -f O(u^) 
where u ~ x — xq. It is assumed that j does not change 
during an AC field cycle. The restoring elastic con- 
stant is then a = d^V / dx'^\x=xa- For ^ parabolic po- 
tential, a = Qfo = const and does not depend on j. 
For a more general form of the potential we take into 
account a cubic term in V{x) in the vicinity oi xq. A 
Labusch parameter obtained from the Taylor expansion 
of the restoring force F = ^dV/dx around xq can be 
written as a = ao(l — (ixo), where /3 is a dimensionless 
constant. Positive (5 describes pinning potential V{x) 
saturating at large a;o, while negative (3 corresponds to 
a potential whose curvature increases with the increase 
of xq. The saturation of a volume pinning force, cor- 
responding to > 0, was discussed by Brandt on the 
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basis of the results of numerical simulations . Satura- 
tion also holds for weak collective pinning. A saturating 
potential provides a smooth crossover to the flux flow 
regime, where the penetration depth diverges (for zero 
viscosity). Negative (3 could be realized for pinning by 
columnar defects. Assuming |/3a;o| <C 1 (small deviation 
of V{x) from parabolic form) we obtain for the vortex 
displacement due to current j: ~ j / ic which gives 
a = ao{l — Pj/jc)- The linear AC response in the pres- 
ence of a Bean current j is then described by: 



(j) - Ai + 



A^c(O) 

1 - Pj/jc 



(1) 



If pinning is weak and magnetic relaxation is fast, as 
is typical for high- Tc superconductors the term f3j/jc 
is small. Equation (|l|) then reduces to: 



AA = A(j)-A(0) = /3 



A^(0) J 
2A (0) 



(2) 




FIG. 1. Schematic illustration of the effective pinning po- 
tential. Solid and dotted lines depict non-parabolic and 
parabolic potentials, respectively. Filled circles represent vor- 
tex position in different situations described in the text. 

Equations (|l|) and (||) predict that a non-parabolic po- 
tential results in an explicit dependence of X on j. For 
a saturating potential (/3 > 0), A(j) is predicted to be 
larger than X{j = 0). For a potential in which the curva- 
ture increases with x {(3 < 0) X{j) would be smaller than 
A(j = 0). Finally, A should depend upon the ratio j/jc 
rather than on j alone. We now address each of these 
three points experimentally and conclude that in typical 
high-Tf, superconductors, X{j/jc) always increases with 
the increase of a ratio j/jc thus providing evidence for a 
saturating pinning potential. 

The penetration depth was measured with an 11 MHz 
tunnel-diode driven LC resonator [^|^ operating in a '^He 
refrigerator. An external DC magnetic field (0-7 kOe) 
was applied using a compensated superconducting mag- 
net. An AC excitation field (~ 5 mOe) was applied paral- 
lel to the DC field. For platelet superconducting samples 



with both AC and DC magnetic fields applied perpen- 
dicular to largest face, an analysis was given in Ref. [||. 
When A is much less than any sample dimension, the 
resonance frequency shift A/ = /(T) — /(T,„i„) of the 
oscillator is proportional to a change in the penetration 
depth, AA = A(r) - X{T„,in) via A/ = -GAA, where 
G is a calibration constant The sample could be 

moved in and out of the sense coil, thus allowing for both 
in situ background subtraction and calibration Q. The 
noise level, A///o « 10~^/-\/ Hz permitted a resolution 
of AAqc < 0.5 A for typical samples. 

The behavior to be described was observed in a vari- 
ety of superconducting single crystals including YBaCuO 
{Tc w 93 K), BiSrCaCuO {Tc w 90 K), PrCeCuO 
{Tc w 21 K), and (BEDT) organic {Tc « 12 and 9.5 K) as 
well as conventional Nb polycrystalline sample {Tc = 9.27 
K). These properties of the pinning potential persist de- 
spite large differences in Tc, symmetry of the order pa- 
rameter, anisotropy and the sign of charge carriers. 
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FIG. 2. Curve 1 2 ~* 3: Hdc = 260 Oe was applied at 
T — 1.5 K. Then, sample was warmed up and cooled down; 
Curve 4 — » 5: field was applied at T = 12 K and sample was 
cooled down; Curve 5 ^ 4 2 — > 3: after previous step, 
sample was warmed up and finally cooled again to T = 1.5 
K. 

Observation of thermal irreversibility: Figure || 
presents a series of temperature sweeps corresponding to 
different ratios oi j / jc{T) for a BSCCO single crystal. In 
the first experiment (solid line 1 ^ 2 ^ 3), the sample 
was cooled in zero field to 1.5 K, a DC field of 260 Oe was 
applied (after first ramping to -7 kOe to insure complete 
penetration), and the sample was then warmed (1 — > 2) 
and then cooled back to 1.5 K (2 ^ 3). Figure |l| shows 
the corresponding trajectory of the system in its pinning 
potential well. Once the temperature exceeded a value 
Tirr{H), j{t) relaxed rapidly and subsequent cooling and 
warming traces (2 ^ 3 ^ 2) were perfectly reversible. 
This curve corresponds to motion near bottom of the pin- 
ning well where both j and /3 are small and A(T) reflects 
changes in jc{T) or, equivalently, the height of the well. 
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We refer to the (2 ^ 3 — > 2) curve as the "reversible" 
curve. It is identical to that obtained in a direct field- 
cooled experiment at the same field. Tirr{H) represents 
a vortex " annealing" or low-T irreversibility temperature 
for the described experiment. In our case, Tirr(260Oe) 
= 28 K, which is close to the temperature at which sev- 
eral changes in vortex structure are often observed 
This temperature is well below the usual irreversibility 
(and melting) temperature determined from magnetiza- 
tion measurements 

The irreversible effects are quite small, '--^ 40 Hz - be- 
tween marks 1 and 3, Fig. The full frequency shift 
from the lowest temperature up to Tc is about 10"^ times 
larger. The same sample, measured using a Quantum 
Design MPMS SQUID magnetometer, showed a diamag- 
netic Meissner signal of < 8 x 10~^ emu. Thus, hystere- 
sis in X(T,j) is equivalent to a change in magnetization 
< 10^^ emu. To resolve the shape of a hysteresis curve a 
sensitivity better than 10~® emu is needed. The effect is 
even smaller on smaller PCCO samples, where a similar 
estimation gives ~ 10~^ emu for the magnitude of the ef- 
fect. This sensitivity is difficult to achieve in commercial 
AC susceptometers, but measurements of large ceramic 
high- Tc samples should, in principle, be able to resolve 
the irreversibility. 
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FIG. 3. Time-logarithmic relaxation of the AC penetra- 
tion depth after application of a 500 Oe magnetic field at 2 K 
(lower curve) and after FC in 500 Oe (upper curve). 

Logarithmic relcixation of X{j): If Eq. (|l|) is cor- 
rect, A(j) should relax along with j, due to thermally 
activated flux creep. Figure ^ shows the time depen- 
dence of \{T,H,t) - X(T,H,to) at T = 2 K in BSCCO 
single crystal. The sample was cooled in zero field, after 
which Hdc ~ 500 Oe was applied. Since the temperature 
was constant jc{T) did not change but X{t) decreased 
logarithmically with time. The decrease in X{t) was in- 
dependent of the sign of j (i.e., if the field was reduced 
after ramping up to a large value). By contrast, the flat 
curve shows X{t) after field cooling in 500 Oe to the same 
temperature. No noticeable relaxation was observed. 
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FIG. 4. Curve 1: A 260 Oe magnetic field was applied after 
ZFC at 1.5 K. Curve 2: After ZFC, magnetic field was first 
set to -7 kOe an returned back to 260 Oe. Schematics shows 
corresponding DC profiles of vortex density marked 1 and 2. 



Dependence upon j/jc- In the next experiment, il- 
lustrated in Fig. H, the sample was cooled to T = 12 K 
and, after ramping to -7 kOe to ensure a complete flux 
penetration, a 260 Oe field was applied (point 4 in Fig. 
H) The sample was then cooled to 1.5 K (4 ^ 5) and 
warmed (5 4 — > 2) above Tirr{H). Again, upon subse- 
quent cooling and warming (2 ^ 3 ^ 2), A(T) remained 
on the reversible curve and coincided with the reversible 
curve obtained in the first experiment described above. 
We interpret this measurement as follows. When a DC 
magnetic field is turned on at Tzfc = 12 K, j is deter- 
mined by U{j/jc{T)) = TzFclri{At/to), where to is the 
characteristic relaxation time and At is the experimen- 
tal time-window When the temperature decreases, 
jc increases and the activation energy U{j/jc) increases, 
thus prohibiting vortices from thermally activated mo- 
tion and essentially freezing the value of j. If A were 
dependent only on j, it would have remained constant 
from mark 4 to mark 5, Fig. |^. The fact that it de- 
creases implies a parametric change in the pinning po- 
tential: when jc increases, j/jc decreases and by Eq. [|, 
the Campbell length decreases. When the temperature is 
again increased (5 4) j{T) remains frozen at its value 
originally set at Tzfc (mark 4) and A(r) retraces the 
preceding cooling (4 5) curve reversibly. Above Tzfc 
K, the activation energy U{j / jc) becomes low enough for 
j to once again relax. Then, A(T) decreases until contri- 
bution of the reversible curve becomes dominant. 

Tilt contribution to the signal, A44: In the configu- 
ration used here (magnetic field normal to the conducting 
planes of a thin platelet), the magnetic field penetrates 
both from sides and from the top and bottom surfaces. 
In the Meissner state, this was demonstrated by solving 
numerically the London equation and testing the results 
experimentally In the presence of vortices, penetra- 
tion from sides. An, is mediated by compressional waves, 
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whereas penetration from the top and bottom faces is due 
to tilt waves, A44 For thin enough samples, the lat- 
ter contribution must dominate the former as we now 
show experimentally. In Fig. |4[ the sample was zero- 
field cooled to T = 1.5 K and a 260 Oe magnetic field, 
less than a field of full penetration, was applied (mark 
1 in Fig. ^). This produced a partially-penetrated state 
where only a fraction of the sample interior was filled 
with vortices, as illustrated by solid lines in the Bean 
profile shown in Fig. ^ Upon warming, j decreased and 
vortices filled the sample. If most of the signal comes 
from the top and bottom surfaces, the penetration depth 
must increase as vortices fill the sample. Only when vor- 
tices reach the sample center should the signal decrease 
in accordance with Eq. (|l|) . Exactly this behavior was ob- 
served on the curve starting at point 1 in Fig. ^ The 
sample was then cooled from above Tc in zero field to 
1.5 K. H was then ramped down to -7 kOe and then re- 
turned back to -1-260 Oe, producing the fully penetrated 
state shown by the dotted and solid lines in Fig. |. The 
resulting curve beginning at 2 is shown in Fig. ^. A(r) 
decreased monotonically upon warming, indicating that 
tilt motion of vortices penetrating the top and bottom 
surfaces constituted the primary signal contribution. 



1000 



500 



< 




e«»i, 



260 Oe 
260 Oe" 



20 30 40 

T(K) 

FIG. 5. Comparison of AA(r) for flux entry and exit. 
Closed symbols: magnetic field was ramped up from - 7 
kOe to -1-260 Oe (flux entry). Open symbols: field was 
ramped down from +7 kOe to -f260 Oe (flux exit) and the 
sample was warmed-then-cooled. Schematics shows the cor- 
responding profiles of vortex density. 

Although the Campbell length in a uniform state de- 
pends upon i?, the observed hysteresis depends explicitly 
upon j /jc- (A weak implicit _B— dependence may arise 
through j/jc{B)). This fact is demonstrated in Fig. |^ 
where we compare A(T) measured at the same final value 
of applied field, H = 260 Oe applied at 1.5 K. The solid 
symbols correspond to the initial application of a - 7 kOe 
magnetizing field while the open ones correspond to + 7 



kOe magnetizing field. Both fields were then ramped to 
H = +260 Oe before measurements began. These field 
changes were sufficient for full flux penetration. The cor- 
responding flux profiles are shown schematically in Fig. 
||. The value of B throughout the sample was much dif- 
ferent for these two starting conditions, but the thermal 
hysteresis curves are essentially identical, demonstrating 
that the hysteresis in Ac is directly related to j/jc, not 
B. The small diflFerences at elevated temperatures could 
be attributed to the differences in j/jc{B). 

In conclusion, we have demonstrated that thermal hys- 
teresis of the AC penetration depth can be attributed to 
changes in j/jc according to Eq.(0). The effect is observ- 
able only for a non-parabolic pinning potential that, as 
follows from our measurements, saturates with displace- 
ment. We identify a new low temperature irreversibility 
line above which thermal hysteresis vanishes. A com- 
prehensive study of TirriH) in conventional and layered 
superconductors will be presented in a forthcoming pa- 
per. 
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